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ABSTRACT 

It is shown, that mechanical alignment can be efficient for suprathermally rotat- 
ing grains, provided that they drift with supersonic velocities. Such a drift should 
be widely spread due to both Alfvenic waves and ambipolar diffusion. Moreover, if 
suprathermal rotation is caused by grain interaction with a radiative flux, it is shown, 
that mechanical alignment may be present even in the absence of supersonic drift. 
This means that the range of applicability of mechanical alignment is wider that it 
is generally accepted and it can rival the paramagnetic one. We also study the latter 
mechanism and reexamine the interplay between poisoning of active sites and des- 
orption of molecules blocking the access to the active sites of H2 formation in order 
to explain the observed poor alignment of small grains and good alignment of large 
grains. To have a more comprehensive picture of alignment we briefly discuss the 
alignment by radiation fluxes and caused by grain magnetic moments. 
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1 INTRODUCTION 

The polarimetry of the interstellar medium (hereafter ISM) 
testifies that dust grains are aligned. The original explana- 
tion of this phenomenon in Davis & Greenstein (1951) as 



being due to par amag netic relaxation was later criticized in 
Jones & Spitzer ( 1967[ ) as inadequate unless enhanced imag- 
inary part of grain magnetic susceptibility was assumed. At 
the same time, alignment involving purely mechanical pro- 
cesses, pioneered in Gold (1951, 1952), was shown to have 
its own problems (see Davis 1955, Purcell 1969, Purcell & 
Spitzer 1971). A new important contribution to the field was 
done by Roberge & Hanany (1990), who proposed an inter- 
esting idea, that grain alignment can be caused by ambipolar 
diffusion. A comprehensive study of alignment of thermally 
rotating grains through the ambipolar diffusion is given in 
Roberge, Messinger, & Hanany (1995). Such an alignment 
is likely to be important for molecular clouds. However it is 
possible to show, that supersonic ambipolar diffusion cannot 
persist over extensive regions of the ISM. 

A really profound step in elucidating the dynamics of 
the ISM dust grains was done by Purcell (1975, 1979) who 
introduced a concept of supratherma l rota tion. However, 
it was shown in Spitzer & McGlynn (1979) that this only 



marginally improves the efficiency of paramagnetic align- 
ment for high rates of resurfacing. In fact, specialists in the 
field believe that 'a relaxation mechanism orders of mag- 
nitude more efficie nt than that given by normal paramag- 
netism' is required (Chlewicki & Greenberg 1990). 
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This paper was preceded by our earlier papers devoted 
to the same issue. First, in Lazarian (1994) (henceforth Pa- 
per I) we provided a quantitative description of mechanical 
alignment of non-spherical grains. However, at that point it 
was unclear whether grains can be aligned mechanically, if 
they rotate at suprathermal velocities. Therefore, only align- 
ment of thermally rotating grains was discussed. Then, in 
Lazarian (1995a) (henceforth Paper II) we addressed the 
problem of paramagnetic alignment of suprathermally ro- 
tating grains and showed that diffusion of oxygen over grain 
surface and existence of potential barrier for H2 formation 
(Tielens & Allamandola 1987) bring into being the grain 
critical size l cr and the critical number of active sites v cr . 
It was found out that for grains with sizes less than l cr 
suprathermal rotation due to H2 formation is suppressed 
and that the poisoning of active sites is accelerated when 
the number of active sites becomes greater than u cr . How- 
ever, only the case v < v cr was quantitatively discussed and 
no processes of desorption were invoked; this is a disadvan- 
tage of the latter study. 

The aim of this paper is to determine roughly the rela- 
tive importance of processes, that can provide alignment of 
suprathermally rotating grains. Therefore the contribution 
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is rather miscellaneous. First of all, in Sect. 2 we remind 
our reader of the Purcell's concept of suprathermal rotation. 
Then, we address in Sect. 3 mechanical processes, which are 
usually believed to be incapable of aligning suprathermally 
rotating grains. We show that this is not true and discuss a 
process that can provide efficient alignment of such grains. 
In Sect. 4 we address the paramagnetic alignment and dis- 
cuss poisoning of grain active sites when v > v cr and the 
processes of desorption of molecules blocking the access to 
active sites. To provide a more comprehensive picture of 
alignment of suprathermally rotating grains, we touch upon 
the alignment of helical grains, grains subjected to radiative 
flows, as well as alignment due to grain magnetic moments. 



2 SUPRATHERMAL ROTATION OF GRAINS 

Here we discuss the suprathermal rotation due to formation 
of H2 molecules on the ISM grains. The number of molecules 
ejected per second from an individual site is ~ 7iZ 2 nErUiz/ _1 , 
where 71 is the portion of H atoms absorbed by a grain of 
diameter I while v is the number of sites. A recoil from a 
molecule departing the grain is rriHi'UH'i = (2niH2£ii2) 1 ' 2 , 
where Eh2 is t he kinetic energy of H2 molecule (0.2 eV) 
( Williams 1988 ). Then the average squared z component of 
torque is 

2 

([Mz] 2 ) « ^l 6 n 2 H m H2 v 2 1 E H 2iy- 1 . (1) 



Thus an individual grain spins up to high angular ve- 
locities, limited only by friction forces. As a result, the 
attainable value of angular velocity f2 is proportional to 
([M z ] 2 ) 1/2 k-\ where k is the constant of rotational friction 
~ Izt^ . Note, that I z is the z component of the momentum 
of inertia and td is the rotational dumping time. The latter 



td £Z Ctn 



c 
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(2) 



where t m is the time that takes a grain to collide with 
gaseous atoms of the net mass equal to the mass of the 
grain, V and q 3 are grain volume and density respectively; 

5 is the grain cross-section for a supersonic flux, the coeffi- 
cient C = 0.6 is a precise result for a spherical grain (Purcell 

6 Spitzer 1971) and $ is the mass flux which is v\nm a for 
supersonic vi. Therefore assuming 71 = 0.2 and ng = n, 
one obtains fi ~ 2 ■ 10 8 s"\ for v = 10 2 . 

Some grains are not subjected to the torques caused by 
H2 formation. For instance, it was shown in Paper II that 
those are aromatic carbonaceous grains, as H2 molecules are 
being formed is states of low excitation over their surface 
(Duley & Williams 1993). Nevertheless, such grains can ro- 
tate suprathermally, for instance, due to variations of the 
accommodation and/or photoelectric emission coefficients 
(see Purcell 1979). Further on, we will show that unlike 
suprathermal rotation due to H2 formation, some of those 
processes may result in alignment on their own. 

Our estimates of Q above are relevant to the long-lived 
spin-up (Spitzer & McGlynn 1979). For short-lived spin- 
up, namely, the correlation time of Purcell rockets ti, is 
much less than td our estimate of Q must be multiplied 
by yth/td- Speaking about Purcell alignment we will be 
mostly concerned with the long-lived spin-up. Paramagnetic 



alignment corresponding to the short-lived spin-up does not 
differ much from the Davis-Greenstein one corresponding to 
the enhanced temperature of grain rotation KT gas , where 



K 



h([M z ] 2 )t d t L 



(3) 



It was pointed out by Bruce Draine (private communica- 
tion), that such an enhancement can be essential in molec- 
ular clouds where temperatures of gas and grains are very 
close. 

In fact, suprathermal rotation can arise not only from 
H2 formation, but also by the variations of the accomoda- 
tion coefficient. We discuss these such processes in Section 5, 
as we show, that such variations can cause mechanical align- 
ment of its own irrespectively from the action of the para- 
magnetic relaxation. 



3 MECHANICAL ALIGNMENT 
3.1 Supersonic motions 

To produce mechanical alignment grains should drift super- 
sonically in respect to gas. In Paper I it is showed that 
Alfvenic perturbations are efficient in providing supersonic 
drift perpendicular to magnetic field lines. Two main phe- 
nomena may be invoked: Alfvenic waves and ambipolar dif- 
fusion. 

It is generally accepted, that supersonic drift by am- 
bipolar diffusion is present within sufficiently strong shocks 
(see fig. 1 in Pilipp et al 1990). However, this cannot be the 
cause of alignment for the majority of the ISM grains. In- 
deed, it was pointed out by B. Draine, that the dissipation 
due to ion-neutral streaming with velocity Vi n is approxi- 
mately 



G — 0.5nni(aTVT)m a Vi, 



(4) 



which for supersonic Vi n = 1km s _1 provides a dissipation 
rate per H atom G/n ~ 1.7 x 10 _25 erg s^ 1 , which is well in 
excess of the energy that supernova is likely to eject into the 
turbulent motions. In other words, the ambipolar diffusion 
may be important for molecular clouds and localized regions 
in within diffuse clouds, but for the large scale alignment 
other processes must be invoked. 

It was hinted in Paper I that large grains are inertial and 
therefore should decouple from ionized gas in the course of 
high frequency Alfvenic oscillations. Although this process 
of alignment does not entail strong dissipation, it is possible 
to show that the limitations imposed on grain size by both 
damping of high frequency Alfvenic waves and grain charge 
are rather strict. 

Indeed, if a grain without electric charge is placed 
within partially ionized gas subjected to Alfvenic perturba- 
tions, e.g. with velocity v = vo sin u At, it is e asy to s ee that 
the amplitude of grain drift scales as Vo/yl + p 2 , where 
\l » 0.4(p n /p g )(vo/(luiA)) and p n and p g are, respectively, 
gas and grain densities. Therefore [i < 1 provides the first 
constrain. In its turn, Alfvenic frequency cannot be arbi- 
trary, as the waves of frequency higher than that of colli- 
sional frequency for ions are critically damped (see McKee 
et al 1994). The corresponding value of uj cr provides the 
minimal value of fi. Another constrain stems from the grain 



© 0000 RAS, MNRAS 000, 000-000 



Alignment of suprathermally rotating grains 3 



charge, as if lo a is less than the grain Larmour frequency, the 
drift velocity scales as vq(uja/wl)- For a given grain size only 
oscillations within the range ujl < uja < ui cr are efficient in 
the sense of providing supersonic drift. Therefore, even if 
the overall non-thermal line broadening of emission lines is 
observed, it is not clear a priori for a given grain, that this 
inequality is fulfilled and that the Alfvenic motions within 
aforementioned range of frequencies are supersonic. The es- 
timates in Lazarian (1995b) show, that for the ionization 
ratio 10 -4 it is difficult to expect grains with radii less than 
10~ 5 cm to experience supersonic drift in the diffuse ISM, 
if grains are collisionally charged (see Spitzer 1978, Draine 
& Sutin 1987). It is gratifying, that such a dichotomy is ob- 
served and grains of sizes less than 10~ 5 cm are not aligned 
(Kim & Martin 1994, 1995). 

Radiation pressure can also cause grain drift. For this 
drift to be supersonic the following inequality QextPrad > 
4:P g as should be satisfied (Purcell 1969), where Q ex t is the 
ratio of the optical to geometrical cross sections, while P ra d 
and P g a a are, respectively, radiation and gaseous pressures. 
For hydrogen with density 10 cm -3 and temperature 100 K, 
Pgas ~ 10~ 13 dyn/cm -2 . Usually Q ex t < 1 and this means 
that it is difficult to obtain sufficient radiation pressure to 
account for large-scale ISM polarization. Thus supersonic 
drift due to radiation pressure is confined to regions in the 
vicinity of bright sources. As grains carry charge, their mo- 
tion under the radiation pressure is usually constrained to 
following magnetic field lines (Spitzer 1978, p. 202). It is 
shown in Habing et al. (1994) that grains should move with 
supersonic relative velocities within the outflows around cool 
giants (carbon-stars, Mira variables and OH/IR stars) due 
to radiation pressure. The velocities of metallic grains should 
be greater than those of dielectric grains in such circum- 
stances (Il'in 1994) and this may have observational conse- 
quences in view of mechanical alignment. 

To summarize, supersonic drift is an essential feature of 
interstellar grain dynamics. In the vicinity of stars this drift 
takes place mostly due to radiation pressure and happens 
along magnetic field lines. Within molecular and atomic 
clouds at sufficient distances from bright sources, this drift 
takes place due to Alfvenic perturbations and grains move 
mostly perpendicular to magnetic field lines. 



3.2 Cross section alignment 

Grain alignment can be quantitatively described by the fol- 
lowing measure: 

OJ = \ (cos 2 / 3 " |) . (5) 

where (5 is the angle between the angular momentum and the 
alignment axis. The latter can coincide with the direction of 
a corpuscular flux, if precession on the time scale of align- 
ment is negligible, or with the direction of magnetic field, if 
the period of grain precession is much smaller than the said 
time scale. For suprathermally rotating grains the axis of 
major inertia is tightly coupled to the angular momentum 
(Purcell 1979). Therefore for oblate spheroidal grains, crj co- 
incides with the Rayleigh reduction factor a introduced in 
Greenberg (1968). For prolate spheroids it is easy to show 
that a — — 0.5a j. 



One of the major questions that we need to answer here 
is whether suprathermally rotating grains can be aligned by 
a supersonic flow. A naive answer would be "no" as such 
rapidly rotating grains resemble gyroscopes and therefore 
should not be sensitive to gaseous bombardment. However, 
this answer ignores the role of crossovers. Further on we 
will show, that when both random and regular torques are 
dominated by recoils caused by Kb formation, an efficient 
alignment of non-spherical grains is possible. 

We remind our reader, that during crossovers grain 
angular velocity becomes close to zero and a substantial 
randomization of J is present. A comprehensive theory of 
crossovers is given in Spitzer & McGlynn (1979) and further 
on we will use the results of this study. 

Let us assume, for simplicity, that randomization is 
complete during a crossover, which is apparently true for 
sufficiently small grains (see Paper II for an explicit expres- 
sion of the disorientation parameter). An incomplete ran- 
domization only alters the time of alignment. For instance, 
grain may need not one, but several crossovers to come to 
the state uncorrelated with an initial one. This, however, 
does not change anything for our further treatment. 

If th is a characteristic correlation time for the existence 
of the sites of H2 formation, the mean time back to crossover 
is given by 

t* » 1.3(t d + t L ) (6) 

(Purcell 1979). Usually, it is assumed that the time of the 
existence of active sites is the time necessary for accreting a 
monolayer of refractory material (Spitzer & McGlynn 1979). 
Another process that can determine ti is poisoning of the 
active sites (see Paper II and Sect. 4 of this paper). However, 
for our simplified treatment we do not distinguish between 
the two processes limiting the life time of Purcell's rockets. 
Indeed, what is important for us is that, for a chosen grain, 
th is the time of accreting of N heavy atoms. The latter 
value depends on particular processes involved, but we will 
see that this is not critical for the mechanism below. 

Consider at first a toy model, namely, assume that the 
axis of suprathermally rotating oblate grain can have only 
two positions, namely perpendicular (_L) and parallel (||) to 
the axis of a gaseous flux. In this model, grain axis stays for 
the time t x in one position and then undergoes the crossover 
and has equal chances to get either the same or the perpen- 
dicular alignment (disorientation is complete!). In the ab- 
sence of the gaseous flux time scales t x \\ and t x ± are the 
same and, naturally, there is no preferential position. How- 
ever, if the gaseous supersonic flux is present, td given by 
Eq. (^) is inversely proportional to gas-grain cross sections 
S n \\ and S n ± in the two positions. Moreover, it is natural 
to assume that the number of heavy atoms adsorbed by the 
grain is proportional to the overall number of atoms strik- 
ing it. Thus tL and t x are both inversely proportional to the 
cross section. Assuming that the time of the crossover is neg- 
ligible as compared with the time of a spin-up, we conclude 
that for an individual grain the time averaged probability 
of finding the grain in a position _L or || is inversely propor- 
tional, respectively, to S n ± and S n \\. To find the constant of 
proportionality one needs to recall, that the probability of 
finding the grain in either of two positions is unity. Accord- 
ing to the ergodic hypothesis this probability coincides with 
the ensemble averaged one. 
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The above considerations can be generalized. It is easy 
to see, that for a continues distribution of axis positions the 
probability of finding grain axis within at a particular angle 
is inversely proportional to the cross section corresponding 
to this angle. 

It was assumed above, that both tL and td are controlled 
by the same process, namely, by the interaction of a gaseous 
flow with a grain. One may imagine situations when tL and 
td are controlled by different physical processes. If one of the 
times tL or td is much less than the other, ^ the longest of 
the two controls the alignment. For instance, for tL C ^ we 
deal with alignment by friction, which may remind one the 
the idea suggested in Salpeter & Wickramasinghe (1969). 
However, these processes are different, as we deal here with 
suprathermally rotating grains, while Salpeter & Wickra- 
masinghe discussed alignment of grains with enhanced rota- 
tional temperature. It is possible to show, that the alignment 
in the latter case is rather marginal. In the opposite limit- 
ing case, namely, ti ^> td we show in Sect. 5 that alignment 
due to photodesorption may be efficient. We do not dwell 
upon all these interesting possibilities here, as believe, that 
the alignment of suprathermally rotating grains under the 
simultaneous action of several processes deserves a separate 
study. 

Consider the alignment of grains subjected to Alfvenic 
perturbations. If, for the sake of simplicity, we approximate 
grains by thin discs the cross section will vary as 



S„ — Tvr sin ipcos tp\ 



(7) 



where ip is the angle between the disc axis and magnetic 
field, V is an angle in the plane perpendicular to the mag- 
netic field and |..| denote that we take the absolute value 
of the trigonometric functions. To account for the Larmour 
precession, we should perform averaging over ip. Therefore, 



(cos 2 ip) = C 



2k 



dip 



cos tp dtp 



(8) 



io \ cos ^\j 
where the normalization constant C is 

(ri^r* (9) 

Jo \cost/j\ 

In short, it is easy to see that (cos 2 <p) = 0.5 and the Rayleigh 
reduction factor (Greenberg, 1968), is equal to 0.25. Grain 
long axis tends to be perpendicular to magnetic field, but 
the alignment is not perfect (compare Paper I). 

In another important case, when flakes stream along 
magnetic field lines, similar computations provide a — —0.5, 
which corresponds to the perfect alignment with grain 
long axis along magnetic field. More efficient alignment for 
streaming along magnetic field as compared with Alfvenic 
perturbations is a consequence of the fact, that in the for- 
mer case the Larmour precession does not change the grain 
- gas cross section. 

It easy to see, that this type of alignment is efficient 
for flakes, as the cross section difference averaged over the 
period of rotation is maximal for such grains. Contrary to 
this, the difference in cross sections is not large for prolate 
grains. For instance, for needles the ratio of the maximal to 

t We should bear in mind that Eq. (JsJ) is an approximate one and 
it is checked in the range of 0.1 < t L /t d < 10 (Purcell 1979). 



minimal averaged cross sections is just 7r/2 and the a for the 
most favorable conditions (a flow parallel to the magnetic 
field lines) is about 0.07. Therefore the alignment of prolate 
grains is marginal due to the mechanism. 

It was implicitly assumed above that the angular mo- 
mentum associated with gaseous bombardment is not im- 
portant during a crossover. It is possible to show that, this 
is true unless grain drift velocity is comparable with ve- 
locity of H2 molecules and/or accommodation coefficient is 
substantially different from unity and/or atomic hydrogen is 
largely converted into molecular form. A detailed discussion 
of the effects of gaseous bombardment is given in Lazarian 
(1995c). 



4 PURCELL ALIGNMENT 

We will call paramagnetic alignment of suprathermally ro- 
tating grains the Purcell alignment to distinguish it from 
the Davis-Greenstein mechanism acting on thermally rotat- 
ing grains. Various aspects of the Purcell alignment were 
addressed in Paper II. Here we extend one aspect of the 
aforementioned analysis, namely, the one dealing with poi- 
soning of active sites. 

Back in Paper II it was established that a critical num- 
ber of active sites exists. This number, v cr , is the mean num- 
ber of sites with chemically adsorbed H atoms, that a hy- 
drogen atom arriving to the grain surface can visit before 
reacting at any of these sites. It is easy to see that, if the 
number of active sites is less than v cr , there cannot be on 
average more than one active site, which is left empty since a 
recent H2 formation. However, if the number of active sites 
v is greater than u cr , the number of empty sites scales as 
u/v cr for v cr 3> 1. These empty sites are the primary tar- 
gets of oxygen atoms hopping over grain surface and it was 
observed in Paper II that for v > v cr poisoning increases. 
Therefore only the case of v < v cr was discussed. Here we 
study the case of v > v cr . Similarly to Paper II we assume 
that oxygen is being immobilized on hydrogenation (Leitch- 
Devlin & Williams 1984, Williams, private communication). 

If vjv cr is the expected number of empty active sites, 
the probability of an oxygen atom to fill any of them as 
a result of an individual hop is v/{v cr N p h), where N p h is 
the number of sites of physical adsorption. Therefore the 
probability of an empty site not to be filled in one hop of an 
oxygen atom is (1 — v / '(y cr Nph))', the same probability for 
mi hops is 



exp 



v rn h 

Vc r N p h 



(10) 



Therefore the characteristic time of poisoning of v/2 
active sites is [[] 

1 



t p — vto 



(11) 



t We assume that not more than one oxygen atom hop over 
grain surface at any particular time. This may not be true for 
dense cores of molecular clouds, where the concentration of atomic 
hydrogen is low. 
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where to is the time of an oxygen atom arrival at grain sur- 
ce. For -± 
Eq ©' 



face. For -^—-^ lL - *C 1 it is possible to expand the exponent 



N, 



ph 



1 



1 m h v 

2 N ph v c 



O 



( ml 


^1 




V% ) 



(12) 



The number of hops is the ratio of the time required 
to hydrogenate oxygen, which we denote t cn to the time of 
an individual hop tu ~ 1CP 12 exp(Eh/(kT s )) s, where Eh 
is the energy of potential barrier and T s is the grain tem- 
perature. When the number of active sites exceeds u cr any 
hydrogen atom scans only the fraction v cr /u of the entire 
surface before either reacting with another H atom at the 
active site or being trapped by an empty active site. The 
probability, that O atom adsorbed by the grain is present 
over this part of the surface is proportional to v cr /v. There- 
fore the time hydrogenation of oxygen is of the order of the 
v jv cr over the timescale of hydrogen arrival. Therefore 



Vcr to , . T 

Vcr " thiVph 

V t H 



1 



v 2 t H 



2 V cr thNph 



O 



,(13) 



where the ratio to/tn = (jinnvo) / (nova); n H & v H and 
no & vo are the densities & velocities of, respectively, hy- 
drogen and oxygen. Note, that for thermal motions vo/vh 
is equal to 0.25. 

If the ratio n — u/N p h stays constant for different 

grains, it is evident from Eq.(13h that t p at the first ap- 
2 — 

proximation is ~ 0.25^-th^-, i.e. it is independent of the 
number of active sites for v > v cr . This is a new result, 
which was not foreseen in Paper II. This result means, for 
instance, that the estimates of t p obtained there for v = p cr 
are applicable to a wide range of grains with v > v cr . 

Poisoning of active sites is essential for the Purcell align- 
ment. A parameter that enters this theory is the ratio of t x 
given by Eq.(^|) to the time scale of paramagnetic relaxation 
tmag (see eq. (58) in Purcell 1979, also see Roberge et al 
1993). To provide sufficient alignment this ratio should be 
greater than unity (see fig. 2 in Purcell 1979). To obtain this 
for standard values of the ISM parameters |^] one needs to 
assume long-lived spin-up, i.e. tL 3> td (Spitzer & McG- 
lynn 1979) and therefore t x ~ ii- It is easy to check that 
tmag scales as I 2 , where / is a grain size and therefore on its 
own paramagnetic alignment favors small grains (see John- 
son 1982). For v > v cr in the first approximation we have 
obtained that t p scales as 1° (see Eq. (12)). 

Note, that for v < v cr our arguments above are not 
applicable. The corresponding study in Paper II showed, 
that in this case not more than one active site is expected 
to be empty and the time of poisoning is proportional to 
vtoN p h/m,h for mu <S N p h- The number of hops for v < v ar 
is the ratio of time scale of the arrival of a hydrogen atom 
tu and the time scale of the hop of the oxygen atom th- 
Therefore both to and mj, scale as l~ 2 , while both N p u and 
v scale as I 2 . As a result the time of poisoning scales as 
I 4 , which makes long-lived spin-up more probable for large 
grains. 



3 It is argued in Paper II that these standard values may be 
misleading for particular regions of diffuse clouds, but we avoid 
discussing this issue here. 



In general, the life time of PurcelPs rockets tL is the 
minimal of the following time scales: t p and the time scale 
of accreting one monolayer of refractory material. In diffuse 
clouds, as a rule, accreting of ice mantles is suppressed (see 
Tanaka et al 1990). It is also believed, that hydrogenated ni- 
trogen and carbon within diffuse clouds do not form mantles 
either. As the abundance of heavier elements is negligible, 
it is natural to assume that in diffuse clouds tL is controlled 
by poisoning of active sites. 

Above we disregarded photodesorption. The character- 
istic time for this process t pa scales as 1° for smooth grains. 
The situation t p d S> t p corresponds to that studied in Paper 
II. If t p d < t p , photodesorption removes molecules blocking 
the access to active sites quicker that these sites are being 
poisoned. Thus a long-lived spin-up is called into being. 

As a result the following qualitative picture emerges. 
For small grains poisoning dominates, i.e. t p d > t p , and 
therefore the spin-up is short-lived. The time of poisoning for 
v < v cr grows with the size as I 4 and for some size becomes 
equal to t p h- Starting from this size we deal with long-lived 
spin-up; tL grows with the size until v = v cr . After the crit- 
ical number of active sites is exceeded, tL is stabilized at 
the attained level. In other words, marginal alignment is ex- 
pected for small grains and good alignment for large grains; 
this corresponds to observations (see Kim & Martin 1994, 
Kim & Martin 1995). 

Using both results obtained above and in Paper II we 
may attempt to write equations for the dinamics of the num- 
ber of active sites. Indeed, while photodesprption, cosmic ray 
bombardment etc. create (or recover) active sites, poisoning 
removes them. The rate of active site creation per unit area 
A may be assumed constant, while the rate of poisoning is 
different for v < v cr and v > v cr . 

If v < v cr the results obtained in Paper II indicate, that 
the rate of poisoning can be approximated by CN~^ , where 
C — 10 12 exp(—Eh/kT s ) s~ , while this rate is of the order 
Therefore for v < v cr the number 



y ph ■ 



of the active sites changes as 
dt N ph 

while for v > v cr the following equation is valid: 



At 



C- 



ph. 



(14) 



(15) 



Evidentely, Eq. ( |14| ) envisages a linear change of the num- 
ber of active sites. If the rate of desorption is greater than 
poisoning, the number of active sites will increase linearly 
with time, unless all the possible active sites are invoked over 
grain surface, provided, that this number is less than u cr ) 
or alternatively the number of active sites is reached u cr . In 
the opposite case when the rate of desorption is less than 
the rate of poisoning the number of active sites decreases 
with time untill all the active sites disappear. If v > v cr the 
solution of the Eq. ( [l5| ) is as follows: 



-TT7— 1 + Bexp( — -,/#^t) 

ANp h ^ "cr-V N ph I 



(16) 



where B is related to the number of active sites Vi n at t = 
in the following way: 
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C 



•An + 



(17) 



Eq. (|L6|) testifies, that the number of active sites tend to 
stabilize at the level v cr ^J AN p h/ C . Therefore for A a bit 
larger than C/N p h, the number of active sites is expected 
to be of the order v cr in correspondence with a qualitative 
conclusion reached in Paper II. 

Unfortunately, in our ignorance of many parameters in- 
volved, it is extremely difficult to quantify this otherwise 
luring picture. First of all, we have rather uncertain knowl- 
edge of the density of active sites as well as of other pa- 
rameters of grain surface. Then, photodesorption presents 
another problem, as the rates obtained in some laboratory 
experiments (see Bourdon, Prince & Duley 1982) are very 
low. We are not sure either whether the photodesorption is 
driven by UV quanta only or also by 3 /jm radiation as it is 
claimed in Williams et al (1992). 

In view of this ambiguities, we need to treat the picture 
above as a conjecture only. We believe, that further research 
in the field will test this conjecture. 

An interesting feature of Eq ([l3]) that it shows that t p is 
proportional to the number of sites of physical adsorption. 
Therefore fractal grains with large surface area should cor- 
respond to larger N p h- The alignment of fractal grains was 
studied in Lazarian (1995d), where the excess of the physi- 
cal area for such grains over that for smooth grains was in- 
voked in order to improve the efficiency of alignment. This 
study indicated a possibility of a substantial improvement 
of alignment, provided that molecules blocking the access to 
active sites had sufficient mobility. This assumption seems 
to be a shortcoming of the latter study. However, it is pos- 
sible to show that the dependences for alignment measure 
on the fractal dimension obtained in Lazarian (1995d) are 
valid without this assumption just as the result of the de- 
pendence of t p on Nph- Another worry in Lazarian (1995d) 
was, that for fractal grains the majority of H2 formation 
events may take place within narrow pores, and this may 
suppress suprathermal rotation. Our arguments above show, 
that the sites within grain pores are likely to be poisoned 
and therefore H2 formation should take place mainly over 
open surfaces of grains, which are kept clean e.g. due to 
photodesorption. At the same time, surfaces within pores 
provide oxygen with more space for its random walk hop- 
ping. 

To summarize, we have shown, that the Purcell (1979) 
mechanism analogously to the Mathis (1986) one favors large 
grains. To tell these two mechanisms it is advantageous to 
study dependences of degree of alignment on grain tempera- 
tures. Our discussion above indicates, that t p is proportional 
to th, which varies exponentially with temperature. There- 
fore the Purcell alignment, unlike the Mathis one should be 
very sensitive to changes of grain temperature. 

Above we assumed that initially the concentration of 
active sites over grain surface is high. Although this corre- 
sponds to the modern picture of grain chemistry (see Tielens 
& Allamandola 1987, Buch & Zhang 1991), an alternative 
approach to the problem is also possible. Indeed, if the sur- 
face density of active sites may be less than 10 -5 cm 2 and 
than only grains larger than 10 -5 cm are likely to have at 



least one active site (Lazarian 1994b) .p] This idea requires 
further study, but here we would like to point out, that the 
case of high and low density of active sites can be distin- 
guished by temperature dependence. It is easy to see, that if 
the density of active sites is low, while nascent H2O molecules 
are ejected on formation, the temperature dependence of the 
Purcell alignment is suppressed. 



5 OTHER POSSIBILITIES OF ALIGNMENT 

Alignment of helical grains. Grains studied above were 
symmetric. However, real grains may have helicity. In fact, 
grains in Purcell (1975, 1979), that rotate suprathermally 
due to variations of the accommodation coefficient or photo- 
electric emission coefficient are helical. What was omitted in 
the above study is that the two aforementioned effects can 
produce alignment even without paramagnetic relaxation. 
For instance, radiation is, as a rule, anisotropic and this may 
influence the rotation caused by photoelectric emission. 

There is a substantial difference when a helical grain is 
subjected to a isotropic bombardment of atoms or photons 
and when it is subjected to a flux. In the former case, there 
is no difference between H2 formation and other causes of 
suprathermal rotation; with high degree of accuracy it is 
possible to assume, that the regular torque acts only along 
the axis of major inertia (see Spitzer & McGlynn 1979). 
If, however, a helical grain is subjected to a radiative or 
corpuscular flux, the component of torque perpendicular to 
the axis of major inertia can cause precession, which may 
significantly alter the alignment. 

An interesting example of helical grains was discussed 
in Dolginov & Mytrophanov (1976). There, grains subjected 
to regular torques due to the difference in scattering of left- 
and right- circular polarized photons were considered. It was 
found, that if grains, either due to their chemical compo- 
sition or due to their shape, have different cross-sections 
for left- and right-hand polarized quanta (see Dolginov & 
Silantev 1976) their scattering of unpolarized light results 
in the spin-up. The efficiency of this process is maximal for 
wavelengths ~ I. For twisted prolate grains, increments of 



grain angular momentum can be shown (Dolginov & Mytro- 



phanov 1976) to be of the order of 0.25& p hl (n r — 1) , where 
<I>p is the flux of photons and n r is the refraction coefficient 
(~ 1.3). Then the characteristic angular velocity is 



Q : 



1 $ p (n. 
8 nm a v a l 2 



(18) 



which is of the order 10 7 s _1 for the typical concentration 
of photons in the ISM ~ 3 ■ 10 9 cm -2 s _1 and may be- 
come much higher in the vicinity of bright sources. For op- 
timal shape suggested in Dolginov & Mytrophanov (1976) 
the scattering efficiency is sufficient to account for the align- 
ment over vast regions, but one may expect the grains to 
have relatively small deviations towards the optimal form. 
Such grains are believed to be a natural product of evolution 
(Mathis 1990). However, numerical studies of the interaction 
of irregular grains with radiation by Bruce Draine (private 

" I suggested this possibility, but did not treat it seriously untill 
it got support from John Mathis. 
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communication) showed, that grains that do not resemble 
helicies can exibit high efficiency of spin-up due to scatter- 
ing of radiation. These studies can shed the light on the 
problem, whether grain spin-up arising from light scatter- 
ing is an exeption or a rule. Note, that such a spin-up is 
expected to be long-lived one. Indeed, it is not the surface, 
but grain volume, that should be altered substantially to 
cause a crossover event. Such a change is expected to take 
place over timescale much longer as compared with the time 
of gaseous damping, which means that the spin-up should 
be long-lived. For such a spin-up the alignment is nearly 
perfect just due to paramagnetic relaxation. 

However, helical grains can be aligned mechanically on 
the timescale much shorter than the one relevant to the para- 
magnetic alignment. Mechanical alignment of helical grains 
is a complex problem and we are going to subject this issue 
to scrutiny in our next paper. Here we just refer to the study 
in Dolginov & Mytrophanov (1976) where it was shown, that 
grains asymptotically tend to align perfectly with their he- 
licity axes along the direction of magnetic field, if the time of 
precession is much less than that of alignment, and along the 
flux direction, if the opposite is true. In their paper Dolginov 
& Mytrophanov considered a tilted oblate grain which he- 
licity axis coincided with the major inertia axis and a tilted 
prolate grain where the two axis were perpendicular. The 
conclusion reached in Dolginov & Mytrophanov (1976) was 
that the two species subjected to a flux should be aligned 
orthogonally. However, the above study omits the influence 
of internal relaxation. Due to this effect a prolate grain ini- 
tially rotating about its axis of minimal inertia will in a short 
period of time turn to rotate about the axis of major inertia. 
In other words, long axis of prolate and oblate grains should 
be aligned in the same direction. 

To obtain quantitative results applicable to the ISM 
and circumstellar regions one needs to estimate the relative 
importance of this particular type of suprathermal rotation. 
Estimates in Purcell (1979) show that, as a rule, angular ve- 
locity of rotation due to H2 formation should dominate for 
the ISM. Therefore, it seems unlikely that alignment due to 
grain helicity is more important in diffuse clouds in compar- 
ison with the one discussed in Sect. 3. However, the align- 
ment of helical grains may be essential in the vicinity of 
bright sources. 

Alignment due to radiation fluxes. Above we dis- 
cussed both mechanical and paramagnetic alignment of 
suprathermally rotating grains. One may wonder whether 
alignment can be caused by absorption of quanta. Indeed, 
due to absorption of an individual quantum a grain gains h 
angular momentum. This physical process is invoked in Har- 
wit (1970). However, it was shown later in Purcell & Spitzer 
(1971) that if a grain is being subjected to a radiation flux, 
the drift produced by the radiation pressure entails Gold 
alignment, which is far more efficient than the Harwit one. 
It is not difficult to show, that the ratio of squared incre- 
ments of angular momentum for the Harwit process ((A J|r)) 
to the one for the Gold process {(AJq)) is negligible even for 
rather special conditions of alignment discussed in Aitken et 
al (1985). Note, that in the aforementioned study the Gold 
alignment was neglected and therefore a conclusion, that 
Harwit mechanism is efficient was reached. We claim, that 
this inference can be true only, if the Gold alignment is sup- 



pressed. The latter happens, for instance, if charged grains 
are trapped by close loops of magnetic field and cannot be 
accelerated under radiation pressure. In any case, such an 
inefficient process as a Harwit one, is unlikely to influence 
the dynamics of suprathermally rotating grains. 

This, however, does not mean, that radiation cannot 
influence alignment of suprathermally rotating grains in a 
way other than through differential scattering of circular po- 
larized quanta. For instance, we may suggest a mechanism 
based on photodesorption. Indeed, if suprathermal rotation 
is due to H2 formation, we may refer to the toy model dis- 
cussed in Sect. 3, but use a radiation flux instead of gaseous 
one. If this radiation flux desorbs molecules blocking the 
access to active sites, the life time of Purcell's rockets and 
therefore £l will be different for _L and || orientation of our 
grain. It is easy to see, that this alignment is most efficient 
if tL S> td- However, in some cases, e.g. for hot and small 
grains discussed in Purcell & Spitzer (1971), the radiation 
dumping may dominate the gaseous one. In this case the 
grain temperature should depend on the grain orientation 
and th S> td is not required for efficient alignment. How- 
ever, a detailed study of these processes is beyond the scope 
of our present paper. 

Alignment due to grain magnetic moments. To 

make our discussion more complete, we need to mention 
the alignment due to grain magnetic moments. One of the 
causes of these moments can be charge (Martin 1971, Draine 
& Salpeter 1979, Draine & Sutin 1987) over rotating grains 
(Rowland effect). The torques that act on a rotating charged 
grain were studied in Davis & Greenstein (1951), along with 
their famous prediction of paramagnetic alignment. How- 
ever, the influence of the ambient gas was omitted in their 
analysis, and thus the conclusion that the 'distribution must 
be completely independent of the field' was made. This was 
quoted in other sources and therefore may still cause confu- 
sion. 

Consider an isolated rotating charged grain in mag- 
netic field. Its rotation creates a magnetic moment Ai = 
Eie(2c)~ 1 a^ J(/) _1 , where ai is the distance from an ele- 
mentary charge over the grain surface to the axis of rotation. 
This causes grain precession in the magnetic field according 
to the Larmour equation 

rl T B 

-=.MxB = . 9 Jx- (19) 

ea 2 B 

where uo g = ^ 3 ' c[ , which is quite close to the grain Lar- 
mour frequency ujl- In the course of this precession the angle 
between Ai and B does not change, i.e. the precession by 
itself cannot change the energy of the system. The calcula- 
tions in Davis & Greenstein (1951) reflect this fact. However, 
the energy changes due to the interaction of grains with sur- 
rounding gas. Indeed, the impact in the direction of the Lar- 
mour precession arising from grain interaction with an atom 
results in the torque, that cause the precession around the 
axis which is perpendicular to the magnetic field. Eq. ( |l9| ) 
shows that the angle between Ai and B increases. If the 
impact is directed against the direction the Larmour preces- 
sion, this angle decreases The collisions of the first type 

II Note, that in reality we assume a uniform distributions of 
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disalign M and B, while the collisions of the second type 
align the two vectors. The alignment effect prevails as the 
impacts against precession are stronger than those in the 
direction of precession. 

In our arguments above we used the fact, that the Lar- 
mour precession influences the direction, but not the mag- 
nitude of the vector J. Therefore the magnitude of M (i.e. 
\M\) is not altered by the Larmour precession. In fact, in 
terms of grain interaction with magnetic field, the grain 
behaves as a magnetic dipole and the alignment due to 
the Rowland effect is similar to the alignment of magnetic 
dipoles. Such an alignment is not complete due to random 
torques acting upon grains.^] These random torques deter- 
mine the effective temperature T and this temperature influ- 
ences the equilibrium distribution of M . If the suprathermal 
rotation is caused by the variations of the accommodation 
coefficient (see Purcell 1979), this temperature will be the 
mean of the grain and gas temperatures, provided that gas- 
grain collisions are inelastic. In the case of the suprathermal 
rotation caused by H2 formation, T will be of the order of 
En2/k, provided that every hydrogen atom that hits grain 
surface leaves it as a part of H2 molecule. 

Indeed, if the spin-up is very long, then grains should 
behave as magnetic moments whose motion is disturbed by 
stochastic torques due to H2 formation.^] These magnetic 
moments are peculiar in a sense, that the randomization 
of their precession requires time much greater than usual 
damping time. However, if the life time of Purcell's rock- 
ets exceeds this time, the treatment in mu228rv2 should be 
applicable. 

In the extreme case of the short spin-up, when all 
torques are essentially stochastic one may introduce the av- 
eraged magnetic moment and study precession of this mo- 
ment in the magnetic field. Our arguments should be appli- 
cable and the damping of the precession should occur over 
the timescale of the order of rotational damping time. 

The equilibrium distribution of grain axes of major in- 
ertia in magnetic field B is proportional to exp(—£/(kT)), 
where £ — —.MB. It is easy to see that 

£ e ^ 2 lotB ujl 
kf ~ 2c kT ~ uj^ 



(20) 



for a thermally rotating grain and £(kT)~ 1 « LOLLOrto^. 2 for 
a suprathermally rotating grain, lot is the frequency of ro- 
tation corresponding to the temperature T, lo t is the fre- 
quency of the suprathermal rotation and lol is the Larmour 
frequency. 

In fact, grain alignment due to the Rowland effect is in 
no way different from the alignment of ferromagnetic grains 
with high permanent magnetization discussed in Spitzer 
& Tukey (1951). The difference is of quantitative nature; 
ferromagnetic grains have much larger magnetic moments 
and therefore much more susceptible for such an alignment. 
However, even with ferromagnetic grains the alignment was 



atomic impacts and subdivide the impacts as having component 
either in the direction of the precession or opposite to it. 
** Thermal fluctuations within grain material can randomize 
alignment as well, but the influence of such fluctuations for 
suprathermally rotating grains is negligible. 

H Note, that we are speaking about averaged J 2 components of 
the torque. 



shown to be not efficient for diffuse clouds P] (Spitzer & 
Tukey 1951). Therefore alignment due to the Rowland ef- 
fect is negligible for the majority of the foreseeable cases. 
This a fortiori true, as paramagnetic grains obtain much 
larger magnetic moments through the Barnett, rather than 
through Rowland effect. Assuming that lot ~ 10 5 s , 
B m 3 x 10~ 6 G, lo t ~ 10 4 lot, we still fall short by approxi- 
mately six orders of magnitude to produce any measurable 
alignment due to this effect. 

To summarize, magnetic moments irrespectively of their 
nature (the Rowland effect, the Barnett effect, permanent 
magnetization) do produce alignment, but this alignment 
is negligible for suprathermally rotating grains in diffuse 
clouds. 



6 CONCLUSIONS 

Several mechanisms of alignment of suprathermally rotating 
grains were discussed and it was shown that 

I. Paramagnetic alignment stays the strongest candi- 
date to account for grain alignment over vast regions of the 
diffuse clouds. Its modification suggested in Purcell (1979) 
depends on subtle processes over grain surfaces and this pro- 
vides an opportunity to test it. 

II. Alignment of grains due to the cross section effect 
that was introduced in this paper can also be important for 
the diffuse ISM. This mechanism, similar to the paramag- 
netic one, tends to align long grain axis perpendicular to 
magnetic field lines, provided that the grain drift is caused 
by Alfvenic perturbations. 

III. Radiation can drive alignment in the vicinity of 
bright sources even in the absence of supersonic drift. Such 
an alignment is more likely to arise from differential pho- 
todesorption, although in some particular cases anisotropic 
radiation field can also provide the alignment of "helical" 
grains. Therefore additional care is required in interpreting 
the corresponding polarimetric data. 

IV. Magnetic moments of grains both arising from the 
Barnett effect and due to their charge produce alignment, 
but the expected degree of alignment is negligible. 
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